Parallel Multigrid Preconditioning of the Conjugate Gradient Method for Systems of Subsurface Hydrology
نویسنده
چکیده
Parallel preconditioners are considered for improving the convergence rate of the conjugate gradient method for solving sparse symmetric positive deenite systems generated by nite element models of subsurface ow. The diiculties of adapting eeective sequential preconditioners to the parallel environment are illustrated by our treatment of incomplete Cholesky preconditioning. These diiculties are avoided with multigrid preconditioning, which can be extended naturally to many processors so that the preconditioner remains global and eeective. The coarse grid correction which deenes the multigrid preconditioner is outlined and its parallel implementation with the distributed nite element data structure is presented, along with some examples of its use as a parallel preconditioner.
منابع مشابه
Preconditioning Techniques for the Bidomain Equations
In this work we discuss parallel preconditioning techniques for the bidomain equations, a non-linear system of partial differential equations which is widely used for describing electrical activity in cardiac tissue. We focus on the solution of the linear system associated with the elliptic part of the bidomain model, since it dominates computation, with the preconditioned conjugate gradient me...
متن کاملUniversity of Colorado at Denver and Health Sciences Center Preconditioned Eigensolver LOBPCG in hypre and PETSc
We present preliminary results of an ongoing project to develop codes of the Locally Optimal Block Preconditioned Conjugate Gradient (LOBPCG) method for symmetric eigenvalue problems for hypre and PETSc software packages. hypre and PETSc provide high quality domain decomposition and multigrid preconditioning for parallel computers. Our LOBPCG implementation for hypre is publicly available in hy...
متن کاملA Parallel Multigrid Preconditioned Conjugate Gradient Algorithm for Groundwater Flow Simulations
This paper discusses the numerical simulation of groundwater ow through heterogeneous porous media The focus is on the performance of a parallel multigrid preconditioner for accelerating convergence of conjugate gradients which is used to compute the pressure head The numerical investigation considers the e ects of boundary conditions coarse grid solver strategy increasing the grid resolution e...
متن کاملNonsymmetric multigrid preconditioning for conjugate gradient methods
We numerically analyze the possibility of turning off post-smoothing (relaxation) in geometric multigrid used as a preconditioner in conjugate gradient linear and eigenvalue solvers for the 3D Laplacian. The geometric Semicoarsening Multigrid (SMG) method is provided by the hypre parallel software package. We solve linear systems using two variants (standard and flexible) of the preconditioned ...
متن کاملNonsymmetric Preconditioning for Conjugate Gradient and Steepest Descent Methods1
We numerically analyze the possibility of turning off postsmoothing (relaxation) in geometric multigrid when used as a preconditioner in conjugate gradient linear and eigenvalue solvers for the 3D Laplacian. The geometric Semicoarsening Multigrid (SMG) method is provided by the hypre parallel software package. We solve linear systems using two variants (standard and flexible) of the preconditio...
متن کامل